Exact bulk correlation functions in one-dimensional nonadditive hard-core mixtures 
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In a recent paper [Phys. Rev. E 76, 031202 (2007)], Schmidt has proposed a Fundamental Mea- 
sure Density Functional Theory for one-dimensional nonadditive hard-rod fluid mixtures and has 
compared its predictions for the bulk structural properties with Monte Carlo simulations. The aim 
of this Brief Report is to recall that the problem admits an exact solution in the bulk, which is 
briefly summarized in a self-contained way. 
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Perhaps the most successful class of density functional 
theories are based on Rosenfeld's Fundamental Measure 
Theory (FMT) In a recent paper @], Schmidt has 
proposed a FMT for the excess free energy of inhomoge- 
neous one-dimensional nonadditive hard-rod fluid mix- 
tures. As a test of the theory, the FMT predictions 
for the pair correlation functions in the bulk region are 
compared with Monte Carlo simulations, a general good 
agreement being found. On the other hand, notwith- 
standing the merits of the FMT constructed in Ref . Q , 
it presents some limitations that become more important 
as the density and/or the nonadditivity increase. For in- 
stance, it yields non-zero values of the pair correlation 
functions inside the core and predicts a spurious demix- 
ing transition. 

It seems to have been overlooked in Ref. [2| the fact 
that the one-dimensional nonadditive hard-rod problem 
admits an exact solution in the bulk. Actually, any one- 
dimensional homogeneous system is exactly solvable, pro- 
vided that every particle interacts only with its nearest 
neighbors |H, 0,11] . The aim of this Brief Report is to fill 
the gap in Ref. [2[ by presenting a brief and self-contained 
summary of the exact solution, particularizing to binary 
nonadditive mixtures, and comparing with the bulk FMT 
predictions for one of the cases considered in Ref. [2|. 

Let us consider an m-component one-dimensional fluid 
mixture with constant (bulk) number densities {pi',i = 
l,...,m} and interaction potentials <pij{ x ) — < t > ij{~ x ') 
acting only on nearest neighbors. Given a particle of 
species i at the origin, the probability that its £th neigh- 
bor belongs to species j and is located at a point between 
x and x + dx is given by p\j\x)dx, what defines the (con- 
ditional) probability density distribution p\ ) (x) . In par- 
ticular, (x) is the nearest-neighbor distribution. The 
distributions pj^ (x) verify the normalization condition 



and obey the recurrence relation 



PvW = E / dx'pg^ix^ix - x>). 
k=i Jo 



(2) 



Its solution in Laplace space is 



PW(«)=[PW(«) 



(3) 



where (s) is the mxm matrix whose elements (s) 

are the the Laplace transforms of pr9 (x) . 

The total probability density of finding a particle of 
species j, given that a particle of species i is at the origin, 
is obtained as 



Pj9ij( x ) =PiA x ) 



1=1 



(4) 



where gij{x) is the pair correlation function. In Laplace 
space, 



P( s ) = P«( s ) . [|- pW( 



G y (s) = —Pij{t 
Pi 

(5) 

where use has been made of Eq. ([3]). Therefore, the 
knowledge of the nearest-neighbor distributions {p\^ {%)} 
suffices to obtain the pair correlation functions {gij{x)}. 
Note that the Fourier transform hij(k) of the total cor- 
relation function hij(x) = g%j{x) — 1 is simply related 
to the Laplace transform Gij(s) of g%j(x) by hij(k) — 
Gij(ik) + Gij(-ik), where i is the imaginary unit. 

It can be proven that the nearest-neighbor distribution 
possesses the following explicit form 



p^\x)=p 3 K l3 e-^\ 



(6) 



dxpf){x) = 1 



(1) 



where j3 — l/fegT and £ = j3p, fcs, T, and p being the 
Boltzmann constant, the temperature, and the pressure, 
respectively. The Laplace transform of Eq. © is 
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p}p(s)=p j K ij tl ij (s + Z), 



(7) 



where Qij(s) denotes the Laplace transform of e 



2 



To close the problem, one needs to determine the am- 
plitudes Kij = Kji and the damping coefficient £. A con- 
venient way of doing so is by enforcing basic consistency 
conditions. Note first that the normalization condition 
(fT]) for I = 1 is equivalent to 



3=1 



Next, since HnLc^oo gij {x) — 1, one must have 
lim sGij(s) = 1. 



(8) 



(9) 



A subtler consistency condition [4( dictates that 



lim x ->ooP« { x )/p\k \ x ) must be independent of the 



,(i)^ /J 1 ) 

ij 

choice of species i. From Eq. ([6]) this implies that 

K, 



-TTT- = independent of i. 

Kik 



(10) 



Equations (f5]) (p~0|) are sufficient to obtain Ay and £. 
To be more specific, let us consider the case of a binary 
mixture (to = 2). Thus, Eq. §5§ yields 



Gh(s) 



Qn(s) [1- OaaW] + 0? a («) 



piD(s) 



G 22 (s) = 



?22 



(a) [1 



(a)] 



2? 2 00 



Gia(s) 



P2^D(a) 



Qigja) 



(11) 



(12) 



(13) 



where 



Qij (a) = \[(hip)P [ ff (a) = ^pIpJK^ (* + £), (14) 

= [1 - Qu(j)] [1 - 22 (a)] - 0? 2 (a)- (15) 
The behavior of Qij (s) for small s is 

Qij (a) = VpipJ^3 (0 + ^ (0* + ^(a 2 )] , (16) 

where f2L(s) is the first derivative of f2y-(s). Application 
of Eq. 0) yields 



An = 



A, 



1- 0^12^12 (Q 

PM0 ' 



i- Pl K 12 n 12 (0 



P2^ 22 (£) 

Next, Eq. ^ implies 

plKun^io + P 2 2 K 22 n' 22 (0 + 2 PlP2 K 12 n' 12 (0 



(17) 



(18) 



-1. 

(19) 



Finally, Eq. (|T0|) becomes 



K n K 22 =Kf 2 . 



(20) 



Equations ([17)) - pO]) constitute a set of four independent 
equations whose solution gives An, K± 2 , A 22 , and £. 
Inserting Eqs. (JT7D and (JTSJ) into Eqs. (JTOJ) and JUJl one 

gets 



A 12 = 



1 



1 + piLn(£) + p 2 L 22 (^) 



PlP2^12(£) £11 (0 + £22(0 - 2Li 2 (0 ' 



(21) 



1 - pA 12 fi 12 (0 + p lP2 [Q? 2 (£) - O u (OO22 (0] #? a = 0, 

(22) 

where we have called Lij(s) = Q' i j(s)/Clij(s) and p = 
Pi + P2 is the total density. Substitution of Eq. (j2"Tj) 
into Eq. (|22|) yields a single equation for £, which in gen- 
eral is transcendental. Once solved, the coefficients Ku 



are obtained from Eqs. pjp. (fT8|) . and (|21|). The exact 
pair correlation functions are then entirely determined in 
Laplace space through Eqs. (fl"l"|) -- (fl"5)) . 

In the particular case of nonadditive hard rods, one 
has e~^* 3 '( x ) = &(x — <7y-), where 0(x) is Heaviside's 
step function, so that 



Lij(s) 



/(a) = s/PlPjKi 



(23) 



(24) 



The constraint to nearest-neighbor interactions implies 
that (Tij < (Jik+Cjk f° r ah i) fc}< In the binary case this 
amounts to 2<r 12 > max(cr 11 , cr 22 ). The recipe described 
by Eqs. $7}, (JTSJ), JST), and ([221) for the thermo- 

dynamic quantity £ = /3p and the amplitudes A^-, and 
by Eqs. pd |) -(TT5 |) and (f24|) for the structural quantities 
Gij(s) are easy to implement. In order to go back to real 
space and obtain the pair correlation functions (x) one 
can use any of the efficient numerical schemes described 
in Ref. [6(. On the other hand, the simplicity of Eq. (|24]) 
allows one to get a fully analytical representation. Note 
first that 



[0n(a) + 022(a) + 0? 2 (a) - On(a)022(a) 



D(s) 



(25) 

When Eq. (J2SJ) is inserted into Eqs. fTT ]H [T3" ]) . one can 

express Gij (s) as linear combinations of terms of the form 



0n 11 (a)0 22 22 (a)0i 2 12 (a) 



nn+ni 2 /2 



e -a(s+C) 

x (p 2 A 22 )" 22+ " 12/2 , (26) 



where a = rancn + n 22 a 22 + n 12 <r 22 and n = n u + 
ri 22 + Ui2. The inverse Laplace transforms g%j(x) — 
C^ 1 [Gij(s)} are readily evaluated by using the property 



3 -a(s+£) 



(a + f 



-e- ^ X , a) ". 6(s-a). (27) 



(n-1)! 
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FIG. 1: (Color online) Bulk pair correlation functions 
gij(x) for a one- dimensional binary hard-rod mixture with 
0-22/0-11 = 2, 0-12/0-n = 15/8, and pi = p 2 = o 1 " 1 1 /4. The 
solid lines are the exact results and the dashed lines are the 
FMT predictions of Ref. 0. 



terms satisfying 8nn + 16n 22 + 15ni 2 < 80 are needed 
for x < IOcth. Moreover, gij(x) = pj 1 p^\x) — Kije~^ x 
in the first shell, i.e., for cry < a; < o"y + Ay, where 
An = min(crii, 2a 12 - an), A22 = min(o- 22 , 2cti 2 - cr 22 ), 
and A12 = min(<7ii, cr 22 ). 

Let us consider a specific system with C22/011 = 2, 
o'12/o'ii = 15/8, and pi = p 2 = o-f 1 1 /4. The cor- 
responding solution of the transcendental equation for 
£ is £ ~ 2.529640-ii 1 , so that /3p/p ~ 5.05927. The 
numerical values of the amplitudes K%j and the con- 
tact values ffij(cri) are K u ~ 91.5298, #22 ~ 1148.60, 



324.24, 



522 (c. 



22 ' 



7.29382, and 



.912(^^2) ~ 2.82473. The property 011(0^) = 522(^2) is 
common to all the equimolar cases (pi = p 2 ), since then 
Eqs. CE7|) and (TH} imply that K n ^n(0 = K 22 ^22 (£)• 
Figure [T] compares the three exact bulk correlation func- 
tions gi-,( x ) with those predicted by the FMT proposed 
in Ref. 2j. The discrepancies are similar to those found 
in Ref. 2| between Monte Carlo simulations and FMT. 

It must be emphasized that the scheme (I5|)- (TTU|) pro- 
vides the exact bulk correlation functions for a one- 
dimcnsional mixture in the absence of external fields. 
The more general problem addressed in Ref. Q, namely 
the excess free energy as a functional of the inhomoge- 
neous densities, is much more complicated and, to the 
best of my knowledge, its exact solution is not known. 
On the other hand, the exact density profiles pj(x) in- 
duced by external potentials Vj (x) can be obtained under 
certain conditions. The trick consists of assuming that 
one of the species (here labeled as i = 0) has a vanish- 
ing concentration (po = 0) and interacts with the other 
species via the potentials 4>oj (x) — Vj (x). The knowledge 
of the bulk correlation functions gij (x) (with pj — > p^ ulk ) 
can then be exploited to get Pj(x) = p}j ul]t goj(x). The 
important limitation, however, is that Vj(x) must repre- 
sent the potential exerted by a wall that acts only on its 
nearest particles. 

To conclude, it is expected that the exact solutions for 
one-dimensional homogeneous systems derived elsewhere 
[H El IE] and summarized in this paper can be useful as 
benchmarks to construct, test, and refine approximate 
theories like the FMT of Ref. 0. This would allow one 
to gain some illuminating insight into the subtleties and 
difficulties of the problem of interest, which can be help- 
ful in its extension to the more realistic case of three- 
dimensional systems. 
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It is important to realize that if one is interested in dis- 
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